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We investigate the interplay of electron-phonon (EP) coupling and strong electronic correlations 
in the frame of the two-dimensional (2D) Holstein t-J model (HtJM), focusing on polaronic ordering 
phenomena for the quarter-filled band case. The use of direct Lanczos diagonalization on finite lattices 
allows us to include the effects of quantum phonon fluctuations in the calculation of spin/charge structure 
factors and hole-phonon correlation functions. In the adiabatic strong coupling regime we found evidence 
for "self-localization" of polaronic carriers in a (tt, tt) charge-modulated structure, a type of superlattice 
solidification reminiscent of those observed in the nickel perovskites La2_xSrxNi04_|_y. 



Electron diffraction measurements on the nickel 
oxide analogue of the layered high-Tc 214 cop- 
per oxide compound, La2-xNixCu04, reveal a 
(tt, 7r)-superstructure spot at quarter filling, i.e. 
for X — 0.5, which has been interpreted as sign of 
a truly 2D polaron ordered phase [1] . 

To study charge- and spin-ordering phenom- 
ena in such systems, exhibiting besides a substan- 
tial EP interaction strong Coulomb correlations, 
let us consider the planar t-J model with an ad- 
ditional on-site Holstein hole-phonon coupling: 

'^H-t-.J — '^ph+'^ho-ph+'^t-J ^ (1) 

'^ho-ph = -V^P^^o'^ibl + bi) h,. (2) 

i 

Here Hph and 'Hf._j represent the phonon part 
and standard t-J model, respectively. In (2), Sp 
is the local EP coupling constant, uiq denotes the 
bare phonon frequency, hi = 1 — J2a '^Irr^ia^ 
b] (c|^) are the usual phonon (projected fermion) 
creation operators. Of course the HtJM gives 
only a very simplified description of the complex 
hole transport and low-spin high-spin interac- 
tions in the nickelates [2]. 

At quarter-filling and for strong EP interac- 
tions, Lanczos studies of the static HtJM (frozen- 
phonon approximation) yield strong indications 
of a Peierls distorted ground state [3] . 

To discuss non- adiabatic effects preserving 
the full dynamics and quantum nature of the 
phonon degrees of freedom, we perform here a 



direct exact diagonalization (ED) of the HtJM 
on a ten-site (TV) square lattice with at most 
M phonons using a well-controlled phononic 
Hilbert space truncation procedure [4]. Since 
memory limitations impose severe restrictions on 
this method, we study an effective polaronic t-J 
model H'j{lf_j{Ai, 7, 7, r^) as well, which can be 
derived from (1) by applying the inhomogeneous 
modified variational Lang-Firsov (IMVLF) ap- 
proach outlined in [4] . The -I- 2 variational pa- 
rameters take into account static displacement 
field (Zii), dynamic polaron (7), finite density 
(7) and squeezing (r^) effects. We stress that 
the IMVLF-Lanczos approach correctly repro- 
duces the adiabatic and anti-adiabatic, weak- and 
strong EP coupling limits [5] . 

In the numerical analysis of the HtJM, we 
first consider the case of spinless fermions (total 

= 5*^02;), i.e., the electronic correlations are 
neglected. Increasing the EP coupling at fixed 
phonon frequency {huo = 0.8), the smooth varia- 
tion of the charge structure factor 5c (tt, tt) in the 
nearly free polaron state {Sp < 2) is followed by 
a strong enhancement in a "quasi-localized" po- 
laron state indicating the formation of a charge 
density wave (CDW) [see Fig. 1]. This crossover 
becomes suppressed in the non-adiabatic regime 
{huiQ = 3). As can be seen from the insets, at Sp — 
3, i.e. in the CDW-like phase, a larger number 
of phonons is still required to achieve a satisfac- 
tory convergence of the ED data. Including more 
phonons we expect an even more pronounced in- 




Fig. 1: Charge structure factor Sc{q) vs Sp, where 
ED data with different M are compared with 
IMVLF results. The insets show the deviation 
of Eq{M) from the IMVLF-value and the mean 
number of phonons Nph in the ground state, re- 
spectively. All energies are measured in units of 
t and we adopt J — 0.4 throughout. 

crease of Sc{tt,tt) [cf. the IMVLF curve]. 

To visualize the hole-phonon correlations, 
Fig. 2 displays the variation of Cho-ph{\i — j\) = 
(tf'ol^i&j&jl^'o) as a function of Sp for both the 
spin-1/2 and spinless fermion cases. Our results 
clearly show the phonon dressing of the holes ac- 
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Fig. 2: Ch o—ph 

(|z - j|) vs ep at all nonequivalent 
distances \Ri — Rj\/a = 0, l,\/2, \/5. Increasing 
M for Ep > 2, we expect an upturn of Cho~ph{0, 2). 



Fig. 3: Magnetic/charge structure factors Sc/s{q) 
and tp,ef f vs Ep for the 2D Ht JM at quarter filling. 

cording to an AB sublattice structure. Due to the 
possible gain of exchange energy J, the density 
oscillations are weakened for spin-1/2 particles. 

Using the IMVLF-Lanczos method to reach 
the strong-coupling adiabatic regime, we notice, 
as Ep increases, a sequence of transitions from 
nearly free to self-trapped polarons [with less 
mobility cx tp^eff = -Ekin(ep, J)/-Ekin(0, J)] so- 
lidifying into a polaronic superlattice and finally 
to charge-separated (CS) states (ep>3). These 
transitions are accompanied by a change of both 
charge and spin structure factors [see Fig. 3]. 
Note that the spin correlations are significantly 
enhanced (weakened) in the CS (CDW) state. 
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